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A. J. Verbiscer and J. Veverka (1990, Icarus 88, 418-428) 
demonstrated that the single scattering properties of terrestrial 
snows are highly forward scattering in contrast to the single 
scattering properties of icy satellite surfaces, which all have a 
backscattering component to their single scattering functions 
as modeled with B. Hapke’s equations (1981, J. Geophys, Res. 
86, 3039-3054; 1984, Icarus 59, 41-59; 1986, Icarus 67, 264- 
280). This study demonstrates to first order that the differences 
between the scattering properties of icy satellite regoliths and 
their terrestrial counterparts are due to differences in particle 
structure and textures. We find that the particles of the regoliths 
of the icy satellites are composed of a complex mixture of grain 
textures and structures. © 1997 Academic Press 


INTRODUCTION 

The bi-directional reflectance variations of the surfaces 
of icy satellites have been modeled using Hapke’s equa- 
tions (Hapke 1981, 1984, 1986) by a variety of researchers 
(Buratti 1984, 1985, 1991, Buratti and Veverka 1983, 1984, 
Buratti et al 1988, 1990, Domingue et al 1991, 1995, Hel- 
fenstein 1986, Helfenstein et al. 1988, Hillier et al. 1991, 
Verbiscer and Veverka 1989, Veverka et al. 1987). Ver- 
biscer and Veverka (1990) used this model to fit terrestrial 
snow reflectance data in a study aimed at comparing the 
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scattering behavior of terrestrial snows and frost to the 
scattering behavior of water ice in the outer Solar System. 
Their results showed that while terrestrial snows and water 
frost are forward scattering, the surface of the icy satellites 
are measured to be backscattering. Verbiscer and Veverka 
(1990) postulated that the source for this difference might 
be more intricate grain textures caused by the low tempera- 
tures and low gravity conditions characteristic of these 
satellite surfaces. Verbiscer et al. (1990) calculated single 
particle scattering functions for aerial and intimate mix- 
tures of forward scattering ice and backscattering Ceres- 
like material. Their goal was to study the possibility that 
the scattering nature of the icy satellites is dominated by 
this backscattering second component which hides the for- 
ward scattering nature of the water ice surface component. 
Verbiscer et al. (1990) found that the ice itself must be 
backscattering to produce the results obtained with 
Hapke’s model. 

A source of criticism and debate over Hapke’s model is 
its tendency to return backscattering asymmetry parame- 
ters when applied to planetary observations and laboratory 
sample measurements (Mishchenko 1994). However, it 
should be noted here that when Verbiscer and Veverka 
(1990) applied this model to terrestrial snow observations 
it returned forward scattering particle scattering functions. 

The source of the scattering difference between terres- 
trial snows and regolith particles of the icy satellites is not 
well understood. The laboratory work of McGuire and 
Hapke (1995) supports the Verbiscer and Veverka (1990) 
initial hypothesis regarding complicated grain textures. 
McGuire and Hapke (1995) studied how shape and the 
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presence of internal scatterers affect the scattering function 
of individual particles. The purpose of their study was 
to examine how a particle that is large compared to the 
wavelength of incident light (a basic assumption in Hapke’s 
model) scatters light dependent on its physical characteris- 
tics. McGuire and Hapke (1995) found that their particles 
with moderate to high densities of internal scatterers were 
backscattering in nature. However, the sizes of the artificial 
particles in their study were on the order of 1 cm in diame- 
ter, which is much larger than most regolith grains. The 
grain sizes of the water ice/frost component on the icy 
Galilean satellites, for instance, have been hypothesized 
to vary from micrometers on Europa to millimeters to 
centimeters on Callisto based on comparisons of their wa- 
ter ice absorption band depths with laboratory measure- 
ments of water ice/frost spectra (Clark 1980, Calvin et al 
1995). While caution needs to be used when applying the 
findings of the McGuire and Hapke (1 995) laboratory stud- 
ies to modeling results of planetary observations, compari- 
sons with their results are still good for first-order determi- 
nations. 

At present there is no evidence to support that the differ- 
ence in scattering behavior between terrestrial snows and 
icy satellite surfaces is due to differences in grain textures. 
We present a first-order test of the particle structure and 
texture hypothesis by comparing the scattering functions 
of terrestrial snows and icy satellite surfaces with the labo- 
ratory results of McGuire and Hapke (1995). McGuire 
and Hapke (1995) determined that at minimum a two- 
parameter scattering function was needed to model their 
single particle scattering measurements, and a double 
Henyey-Greenstein function modeled their results more 
accurately than a two-parameter Legendre polynomial. 
Hartman and Domingue (1997) investigated the effects of 
the choice of single particle scattering function on the 
Hapke model by determining whether a two-parameter 
or a three-parameter double Henyey-Greenstein function 
better modeled the McGuire and Hapke (1995) measure- 
ments and whether the model predicts significantly differ- 
ent results based on their findings with the laboratory mea- 
surements. They concluded that the three-parameter 
double Henyey-Greenstein function did not provide a sub- 
stantial improvement or difference over the two-parameter 
double Henyey-Greenstein function when used in 
Hapke’s model. However, Hartman and Domingue (1997) 
did find that the choice of single particle scattering function 
incorporated into Hapke’s model did influence the values 
of the remaining model parameters when applied to phase 
curves of the lunar near side. Therefore we have fitted the 
terrestrial snow data of Middleton and Mungall (1952) 
using both double Henyey-Greenstein functions in con- 
junction with Hapke’s model. 

Based on Hartman and Domingue’s lunar phase curve 
results, Domingue and Verbiscer (1997) re-analyzed the 


solar phase curves for the icy Galilean satellites using both 
two- and three-parameter double Henyey-Greenstein 
functions. Their goal was to test whether the hemispherical 
dichotomies in surface structure and texture deduced from 
Hapke’s model were real or a product of the single scatter- 
ing function incorporated into the model. Domingue and 
Verbiscer (1997) found that while the absolute values of 
the model parameters were dependent on the choice of 
scattering function, the trends in the hemispherical dichot- 
omies predicted by this model were independent of the 
chosen scattering function and consistent with the exogenic 
processes modifying these satellite surfaces. Since the ter- 
restrial snow data analyzed by Verbiscer and Veverka 
(1990) were modeled using a single term Henyey- 
Greenstein function, as were many of the icy satellites to 
which they compared their results, we have refitted the 
terrestrial snow data for ease of comparison with the 
McGuire and Hapke (1995) laboratory results and with the 
icy Galilean satellite results from Domingue and Verbiscer 
(1997). In addition we also reanalyzed Rhea’s solar phase 
curve from Domingue et al (1995) with a three parameter 
Henyey-Greenstein single particle phase function for fur- 
ther comparisons with icy Galilean satellite and laboratory 
results. Finally, we compare our results to the Hartman 
and Domingue (1997) analysis of the McGuire and Hapke 
(1995) measurements to infer first-order particle structures 
for both the terrestrial snows and the icy satellite rego- 
lith particles. 

DATA 

The terrestrial snow data analyzed by Verbiscer and 
Veverka (1990) were the in situ goinophotometric mea- 
surements of natural snow and frost from Middleton and 
Mungall (1952). These are the same measurements ana- 
lyzed here. Middleton and Mungall measured the bidirec- 
tional reflectance of six types of naturally occurring snows 
at reflection angles of -80° to +80° and incidence angles 
of 0°, 30°, 45°, 60°, and 75°. Housed in their portable gonio- 
photometer were a collimated light source, which could 
be positioned at 0°, 30°, 45°, 60°, and 75° from the surface 
normal, and a photomultiplier tube that rotated on concen- 
tric tracks from -80° to +80°. The light source and photo- 
multiplier were offset in order to prevent obstruction of 
the light source at very small phase angles. As a result of 
this displacement, a phase angle of zero was impossible to 
measure. For more details on the apparatus used for these 
measurements, the reader is referred to Middleton and 
Mungall (1952). Five of the natural snows were chosen by 
Verbiscer and Veverka (1990) for their study (and there- 
fore this project). These are: (1) newly fallen snow, (2) 
surface frost, (3) rain crusted snow, (4) wind blown snow, 
and (5) settling snow. The sixth type of snow observed by 
Middleton and Mungall, “glazed rain crusted snow,” was 
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TABLE I 

Satellite Phase Angle Coverage 


Satellite 

Data 

source 

0.47-ju.m 

Lead 

o 

H is 
3 

0.55-)Ltm 

Lead 

0.55-/xm 

Trail 

Europa 

Telescopic 

0.2- 11.3 

0.2- 11.3 

0.2- 11.3 

0.2- 11.3 


Voyager 

3.8-119.0 

3.0-119.0 

3.7-104.4 

2.9- 94.3 

Ganymede 

Telescopic 

0.6- 11.4 

1.4- 11.2 

0.4- 11.4 

1.4- 11.2 


Voyager 

4.8-120.3 

20.0-124.0 

4.9- 32.2 

20.0- 38.7 

Callisto 

Telescopic 

0.4- 11.1 

0.6- 11.2 

0.4- 11.1 

0.6- 11.2 


Voyager 

8.3- 50.4 

13.9-140.8 

8.3- 51.7 

14.0-107.4 

Rhea 

Telescopic 

0.3- 6.4 

O 

i 

Os 

0.3- 6.4 

0.9- 6.9 


Voyager 

13.1-135.5 

1.8- 16.4 

13.3-135.6 

1.8- 16.6 


Note. Lead corresponds to leading hemisphere, and trail corresponds 
to trailing hemisphere. All angle ranges are in degrees. 


not analyzed since it had a large specular component which 
Hapke’s equations do not model well. 

The planetary phase curves modeled in this study include 
both hemispheres of the icy Galilean satellites and Rhea at 
0.47 and 0.55 /tm. The Galilean satellite data and modeling 
results are taken from Domingue and Verbiscer (1997). 
The Rhea data and two-parameter double Henyey- 
Greenstein modeling results are taken from Domingue et 
al. (1995). The phase curve data for all the icy satellites 
are a combination of telescopic observations and Voyager 
image data. Table I lists the phase angle coverage available 
for each of these satellites. For a detailed description of 
the data and data processing the reader is referred to the 
above references. 

ANALYTICAL METHODS 

The model used in this project is based on radiative 
transfer. An integral part of this model, or any photometric 
model that attempts to describe scattering from a planetary 
regolith, is the incorporation of a single particle scattering 
function. The role of the single particle scattering function 
is to describe how an average particle or grain from the 
planetary surface preferentially scatters light in one direc- 
tion or another. This knowledge is an important tool for 
providing supporting evidence for various geologic pro- 
cesses such as ion bombardment, where irradiation of an 
ice surface by ions in a corrotating magnetic field can an- 
neal lattice defects and cause the scattering function to 
become more forward scattering (Sack et al 1992). For 
this project we chose two forms of the double Henyey- 
Greenstein function: a two-parameter double Henyey- 
Greenstein function (2P-HG) which is given by 

P 2 (a) = ((1 - c 2 )(l - bl)![ 1 + 2 b 2 cos (a) + bjf' 2 ) 

+ (c 2 (l - bl)l[ 1 - 2 b 2 cos (a) + b |] 3/2 ), 

where a is the phase angle and b 2 and c 2 are the single 


particle scattering parameters, and a three-parameter dou- 
ble Henyey-Greenstein function (3P-HG) whose form is 
given by 

P 3 (a) = ((1 - c 3 )(l - bM 1 + 2 b 3 cos (a) + bjf 12 ) 

+ (c 3 (l - d 2 3 )/[ 1 - 2 d 3 cos (a) + dl] V7 ), 

where a is the phase angle and b 3 , c 3 , and d 3 are the 
single particle scattering parameters. Note that the 3P-HG 
function defaults to the 2P-HG function for cases where 
d 3 = The parameter c in both forms of the double 
HG-function is a measure of the relative amplitudes of 
the forward and backscattering components of the single 
scattering function. The parameters b and d are the widths 
of the forward and backscattering lobes, respectively, and 
are defined as the mean cosines of the scattering angle. 
Therefore the 2P-HG function assumes that the widths of 
the lobes are identical and just the amplitudes vary. 

A modified least squares grid search was used to fit 
Hapke’s disk-resolved equations (Hapke 1984, 1986) to 
the five types of terrestrial snow observations used in the 
Verbiscer and Veverka (1990) study. The resulting model 
parameters are listed in Table II with the corresponding 
rms values. The rms values are defined such that 



where N is the number of observational data points, R is 
the calculated reflectance, r is the observed reflectance, 
and the sum is over the number of points. Figures 1 through 
5 plot our theoretical results versus the data sets. 

When Verbiscer and Veverka (1990) fit these data sets 
they set the opposition amplitude parameter, bo, equal to 
zero since the snow data show no opposition surge. They 
also set the roughness parameter, 0- bar, to zero since the 
types of snows studied here had no visible surface 
roughness. We chose to let 0-bar vary; however, our results 
show that indeed these snows are very smooth, consistent 
with the Verbiscer and Veverka (1990) assumptions. How- 
ever, we also chose to allow both opposition parameters 
to vary over our grid search, since the measurements of 
Middleton and Mungall (1952) contain observations at 
i = 0.0° and e < 1°. We found a few percentage differences 
(up to 2% in some cases) in the rms values obtained when 
we allowed bo to vary compared to when we set bo equal 
to zero. In the model analyses where we set bo = 0 the 
resulting values of the remaining parameters were similar 
to those obtained in the model analyses where bo was 
allowed to vary as demonstrated by comparing Tables II 
and III. The most significant differences were in the single 
scattering function parameters for the newly fallen snow 
and snow frost. These differences predict a slightly less 
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TABLE n 

Hapke Parameters for Terrestrial Snow 



Newly 

fallen 

Snow 

frost 

Rain 

crusted 

Wind 

blown 

Settling 


Newly 

fallen 

Snow 

frost 

Rain 

crusted 

Wind 

blown 

Settling 

w 

0.995 

0.992 

0.986 

0.990 

0.992 

w 

0.995 

0.993 

0.986 

0.990 

0.992 

bo 

1.0 

1.0 

1.0 

1.0 

1.0 

bo 

1.0 

1.0 

1.0 

1.0 

0.5 

h 

0.995 

0.995 

0.995 

0.995 

0.995 

h 

0.995 

0.995 

0.995 

0.995 

0.995 

b 2 

0.449 

0.82 

0.624 

0.632 

0.598 

b 3 

0.445 

0.785 

0.624 

0.632 

0.600 


0.004 

0 

0 

0 

0 

Ci 

0.001 

0 

0 

0 

0 







d. 

-0.704 

0 

0 

0 

0 

frbar 

0 

0 

3 

0 

0 

Abar 

0 

7 

3 

0 

0 

rms 

0.040 

0.025 

0.034 

0.063 

0.033 

rms 

0.041 

0.024 

0.034 

0.063 

0.033 

Note. 

Error bars: w 

= ±0.01, 

bo = ±0.05, h 

= ±0.005, b 2 

= c 2 = b 3 

II 

3 - 

ii 

a 

ii 

±0.02, 0-bar 

= ±2. 





Newly Fallen Snow 



Newly Fallen Snow 



Newly Fallen Snow 



observations 

2P-HG Hapke model solution 
3P-HG Hapke model solution 


FIG. 1. Comparison of the newly fallen snow reflectance measurements (open squares) of Middleton and Mungall (1952) with the theoretical 
modeling results using Hapke’s model incorporating a two-parameter double Henyey-Greenstein single particle scattering function (open diamonds) 
and the theoretical modeling results using Hapke’s model incorporating a three-parameter double Henyey-Greenstein single particle scattering 
function (open triangles). 
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2P-HG Hapke model solution 
3P— HG Hapke model solution 


FIG. 2 Same as described in the legend for Fig. 1, for snow frost. 


forward scattering function for the newly fallen snow for 
the bo = 0 case than for the bo = 1 case. However, the 
differences predict a slightly more forward scattering func- 
tion for the snow frost in the bo = 0 case than for that in 
the bo = 1 case. Our discussions pertain to the bo = 1 
case since there are opposition measurements within the 
data. The opposition width parameter, h , was found to be 
uniformly wide (h — 0.995) for all the snow data set results 
where the opposition parameters were allowed to vary. A 
value of h 1 signifies that there is essentially no opposi- 
tion surge in the snow data. 

A modified least-squares grid search was used to fit 
Hapke’s disk-integrated equation (Hapke 1984, 1986) to 
the disk-integrated phase curves of the three icy Galilean 
satellites and Rhea. The icy satellite data consist of two 
phase curves for both leading and trailing hemispheres, 
one based on observations taken at 0.47 /Am and the other 
at 0.55 /air. Each phase curve was modeled twice, once 


using Hapke’s model with the 2P-HG function and a sec- 
ond time using Hapke’s model with the 3P-HG function. 
The resulting model parameters are listed in Table IV with 
the corresponding rms values. The rms values are defined 
such that 


rms = 


2(A/„ - m„) 2 



where N is the number of observational data points, M„ 
is the calculated magnitude from Hapke’s theory, m„ is 
the observed magnitude, and the sum is over the number 
of points. Figures 6 and 7 plot the theoretical models versus 
observations for the satellite Rhea at 0.47 and 0.55 /Am, 
respectively. Similar plots for the icy Galilean satellites are 
shown in Domingue and Verbiscer (1997). 

The error bars listed in Table IV are based on the grid 
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Rain Crusted Snow 



Rain Crusted Snow 
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legend 
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2P-HG Hopke mode! solution 
3P— HG Hapke model solution 


FIG. 3. Same as described in the legend for Fig. 1, for rain crusted snow. 


size of the modified least-squares fitting routine. The de- 
tails of this fitting routine are described in Domingue and 
Verbiscer (1997). However, the phase angle coverage in 
the satellite phase curves is also important in establishing 
the quality of fit for many of the parameters. A detailed 
error analysis of the fits to the icy Galilean satellites is 
given in Domingue and Verbiscer (1997) which is also 
applicable to Rhea. The final error bars for each of the 
model parameters are given in Table V. 

RESULTS 

Terrestrial Snows 

Hapke’s equations model the terrestrial snow observa- 
tions reasonably well, as demonstrated in Figs. 1 through 
5. Comparisons of the results between the 2P-HG and 
the 3P-HG models are surprisingly similar. Both models 
predict high values of the single scattering albedo (w) which 


are within the error bars (±0.002) of each other. In addi- 
tion, both models predict very smooth values for the mac- 
roscopic roughness parameter, 0-bar. In general, values for 
0-bar lie in the range 0-3°, all approximately within the 
error bars (±2°); however, the 3P-HG model for snow 
frost predicts a rougher surface than the 2P-HG model. 

Both models predict the same value for the opposition 
width parameter, h, for all snow types. The high value of 
h is consistent with the lack of an opposition surge in the 
snow data. Hapke’s model (1986) correlates the value of 
h to porosity within the optically active layer of the regolith. 
However, to obtain an accurate estimate of the regolith 
porosity a reasonable measure of h is needed. Due to the 
lack of the detection of an opposition surge in the snow 
data, the value of h obtained in the modeling results cannot 
be accurately correlated to porosity. 

The opposition amplitude parameter, bo, can be corre- 
lated to the opacity of the regolith particles. It is a measure 
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Wind Blown Snow 


Wind Blown Snow 



Wind Blown Snow 
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2P-HG Hapke mode! solution 
3P-HG Hapke model solution 



FIG. 4. Same as described in the legend for Fig. 1, for wind blown snow. 


of the percentage of incident light which is scattered at 
the surface; thus, a value of 1 signifies that all the incident 
light is scattered at the surface and the particle is opaque. 
The value of bo is known to influence the shape of the 
disk-integrated phase curve out to 30° for which the terres- 
trial snow data sets have corresponding measurements. 
However, the value of bo is influenced by the value of h. 
Since there was no detection of an opposition surge in the 
snow observations and h has a large value corresponding 
to the lack of an opposition surge, no interpretation of the 
quantitative values of bo can be made. 

The single scattering functions from both models are 
plotted in Fig. 8. They both predict highly forward scatter- 
ing behavior for terrestrial snow and water frost, consistent 
with the results of Verbiscer and Veverka (1990), who fit 
the snow data with a version of Hapke’s model which 
incorporated a single Henyey-Greenstein function. No sig- 
nificant differences are seen in the single particle scattering 


function predicted by Hapke's model incorporating the 
2P-HG function versus this same model incorporating the 
3P-HG function. Except for the newly fallen snow data, 
the values of c 2 and c 3 are uniformly zero. When c 2 and c 3 
are zero, the 2P-HG and 3P-HG functions default to the 
single Henyey-Greenstein function used by Verbiscer and 
Veverka (1990). Even the c 2 and c 3 values for the newly 
fallen snow are zero within the error bars. Our results also 
agree with the values for the single scattering function 
found by Verbiscer and Veverka (1990). We attribute this 
lack of a backscattering component to the unique scattering 
nature of snow. 

Satellite Surfaces 

Buratti (1991) modeled both hemispheres of Ganymede 
and Callisto using Hapke’s model to look for hemispherical 
variations in their surface photometric properties. The val- 
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Settling Snow 
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Settling Snow 
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FIG. 5. Same as described in the 


ues quoted by Buratti for the surface roughness parameter 
are commensurate with our Ganymede results, but greatly 
differ from our Callisto results. Our values for the single 
scattering albedo agree with Buratti’s values for Ga- 


Settling Snow 



emission angle (degrees) 


Settling Snow 



emission angle (degrees) 


legend 
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2P-HG Hapke model solution 
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legend for Fig. 1, for settling snow. 


nymede's leading and Callisto’s trailing hemispheres; how- 
ever, they disagree with her values for Ganymede's trailing 
and Callisto's leading hemispheres. We attribute these dif- 
ferences to the different single particle scattering functions 


TABLE III 

Hapke Parameters for Terrestrial Snow with bo = 0 



Newly 

fallen 

Snow 

frost 

Rain 

crusted 

Wind 

blown 

Settling 


Newly 

fallen 

Snow 

frost 

Rain 

crusted 

Wind 

blown 

Settling 

w 

0.998 

0.993 

0.988 

0.992 

0.992 

w 

0.998 

0.990 

0.988 

0.992 

0.993 

b 2 

0.503 

0.716 

0.610 

0.617 

0.594 

h 

0.503 

0.695 

0.610 

0.617 

0.594 

Cl 

0 

0 

0 

0 

0 

C 3 

0 

0.122 

0 

0 

0 







di 

0 

-0.161 

0 

0 

0 

0-bar 

0 

6 

3 

0 

0 

0- bar 

0 

0 

0 

0 

0 

rms 

0.048 

0.025 

0.052 

0.077 

0.036 

rms 

0.046 

0.023 

0.052 

0.077 

0.036 


Note. Error bars: w = ±0.01, b 2 = c 2 = £>3 = d 3 = c 3 = ±0.02, 0-bar = ±2. 
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TABLE IV 

Hapke Parameters for the Icy Galilean Satellites and Rhea 


Leading 


Trailing 


Leading 


Trailing 



0.47 jum 

0.55 fim 

0.47 jim 

0.55 jxm 


0.47 fjim 

0.55 jLim 

0.47 fim 

0.55 gm 

w 

0.922 

0.964 

0.897 

Europa 

0.940 

w 

0.934 

0.964 

0.897 

0.930 

bo 

0.50 

0.50 

0.45 

0.50 

bo 

0.49 

0.43 

0.51 

0.521 

h 

0.0016 

0.0016 

0.0016 

0.0016 

h 

0.0015 

0.0016 

0.0016 

0.0016 

b 2 

0.431 

0.429 

0.43 

0.443 

bz 

0.770 

0.726 

0 

0.083 

c 2 

0.921 

0.887 

0.713 

0.609 

cz 

0.780 

0.945 

0.691 

0.784 

0-bar 

10 

10 

10 

10 

ch 

0-bar 

-0.459 

8 

-0.416 

10 

-0.417 

11 

-0.386 

11 

rms 

0.029 

0.029 

0.039 

0.033 

rms 

0.029 

0.029 

0.035 

0.025 

w 

0.830 

0.945 

0.870 

Ganymede 

0.81 

w 

0.830 

0.930 

0.870 

0.810 

bo 

0.62 

0.86 

1.00 

0.23 

bo 

0.63 

0.91 

1.00 

0.230 

h 

0.003 

0.004 

0.074 

0.074 

h 

0.003 

0.003 

0.074 

0.074 

b 2 

0.282 

0.380 

0.039 

0.307 

bz 

0.090 

0.200 

0 

0.350 

c 2 

0.960 

0.427 

0.989 

0.962 

cz 

0.960 

0.450 

0.820 

0.970 

0-bar 

28 

29 

35 

35 

dz 

<9-bar 

-0.280 

28 

-0.380 

28 

-0.046 

35 

-0.305 

35 

rms 

0.067 

0.036 

0.097 

0.047 

rms 

0.067 

0.036 

0.097 

0.047 

w 

0.740 

0.540 

0.470 

Callisto 

0.550 

w 

0.510 

0.605 

0.530 

0.650 

bo 

1.00 

1.00 

0.27 

0.73 

bo 

0.82 

1.00 

0.50 

0.87 

h 

0.031 

0.031 

0.0277 

0.0277 

h 

0.031 

0.031 

0.0277 

0.0277 

b 2 

0.729 

0.132 

0.432 

0.206 

bz 

0 

0.005 

0.787 

0.988 

c 2 

0.024 

0.949 

0.542 

0.958 

Cz 

0.033 

0.015 

0.489 

0.618 

0-bar 

42 

42 

42 

42 

dz 

0-bar 

-0.694 

42 

-0.687 

42 

-0.400 

42 

-0.238 

42 

rms 

0.093 

0.062 

0.213 

0.066 

rms 

0.092 

0.061 

0.195 

0.065 

w 

0.980 

0.986 

0.927 

Rhea 

0.951 

w 

0.995 

0.996 

0.996 

0.924 

bo 

0.52 

0.46 

0.50 

0.44 

bo 

0.45 

0.44 

0.45 

0.57 

h 

0.0071 

0.0071 

0.0074 

0.0074 

h 

0.0071 

0.0071 

0.0074 

0.0074 

b 2 

0.569 

0.575 

0.567 

0.567 

bz 

0.659 

0.414 

0.352 

0.024 

c 2 

0.336 

0.331 

0.335 

0.323 

Cz 

0.157 

0.171 

0.020 

0.964 

0-bar 

15 

15 

16 

16 

dz 

0-bar 

-0.665 

17 

-0.650 

15 

-0.822 

16 

-0.345 

16 

rms 

0.070 

0.032 

0.081 

0.048 

rms 

0.065 

0.028 

0.079 

0.042 

Note . Error bars: vv = 

±0.01, bo = 

±0.01, /7 = ±0.001, &2 

and b 3 = ±0.005, C 2 and c 3 = 

±0.005, dz = 

±0.005, 0-bar = 

±2. 



incorporated into the models. Buratti used a single term 
Henyey-Greenstein function in her analysis; therefore, her 
results are not easily compared with those stated in this 
paper. Helfenstein et al (1996) modeled Callisto’s surface 
using the same model with a 3P-HG single scattering func- 
tion. Their results differ from those stated in this paper; 
however the error bars quoted by Helfenstein et al (1996) 
place our results close to theirs. The major difference be- 
tween their analysis and ours is that they allowed the oppo- 
sition amplitude parameter, bo, to have values greater than 
1 This will affect the values of the remaining model param- 
eters; therefore in the following discussions we will not 


compare our results with theirs in order to maintain inter- 
nal consistency. 

The values for the single scattering albedo between the 
two models are within the error bars for the satellites 
Europa and Ganymede. The only differences in w seen 
for Europa were between the 0.47-/xm leading hemisphere 
results and the 0.55-jam trailing hemisphere results. The 
only difference in w seen for Ganymede was between 0.55- 
jam leading hemisphere results. The two models returned 
significantly different results for both hemispheres of Cal- 
listo at both wavelengths of observation; however, the lead- 
ing hemisphere of this satellite has no large phase angle 
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phose angle (degrees) 




FIG. 6. The solar phase curves of Rhea’s leading (a and b) and trailing (c and d) hemispheres at 0.47 /im, where the diamonds represent the 
disk-integrated observations, the solid line represents the Flapke model solution incorporating a two-parameter HG function, and the dashed line 
represents the Flapke model solution incorporating a three-parameter HG function. 


data to constrain its phase curve at these important viewing 
geometries. The large phase angle data constrain the values 
of w, 0-bar, and the single particle scattering function pa- 
rameters. There were no consistent trends in the Callisto 
results. The two models also returned significantly different 
values of w for Rhea’s trailing hemisphere. The leading 
hemisphere data set values for w were within the error bars; 
however, there are significant differences in w between the 
two models for the trailing hemisphere data sets. This is not 
surprising since like Callisto’s leading hemisphere, Rhea’s 
trailing hemisphere has no large phase angle data to con- 
strain the model parameter values. 

The value of the surface macroscopic roughness parame- 
ter, 0-bar, showed no significant variation with model in 
the results for the icy satellites. The only other Hapke 
parameter for which there were no significant variations 
with model was the opposition width parameter. 

The opposition amplitude parameter, however, does 


show some significant differences between the two models. 
We define a significant difference in bo to be a difference 
of 0,1 which is based on variations in rms as bo varies 
while the other parameters remain constant. Changes in 
rms of 5% are considered significant. The only satellites 
showing a signifcant difference in bo between the two 
model results were Callisto and Rhea. The significant dif- 
ferences in bo values in the Rhea modeling results were 
for the trailing hemisphere at 0.55 /xm only and can be 
attributed to the scatter in the opposition measurements. 
However, the opposition data for Callisto have remarkably 
little scatter; therefore this difference in opposition param- 
eters is not readily understood. These variations in bo with 
model for these satellites are consistent with the terrestrial 
snow results. 

The single particle scattering functions show strong vari- 
ations with choice of model in the satellite solutions. Fig- 
ures 9 through 12 plot the scattering functions for Europa, 
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TABLE VI 

Isotropic vs Anisotropic Multiple Scattering Model Results for Rhea 


Leading at 0.47 jjl m 


Leading at 0,55 /xm 


Isotropic Anisotropic Isotropic Anisotropic 


w 

0.995 

±0.01 

0.9909 

±0.0009 

0.996 

±0.01 

0.9935 

±0.0002 

bo 

0.45 

±0.01 

0.739 

±0.8337 

0.44 

±0.01 

0.549 

±0.2399 

h 

0.007 1 

±0.001 

0.0539 

±0.232 

0.0071 

±0.001 

0.0041 

±0.0029 

b 3 

0.659 

±0.005 

0.4828 

±0.07 

0.414 

±0.005 

0.3394 

±0.1263 

<3 

0.157 

±0.005 

0.17 

±0.09 

0.171 

±0.005 

0.17 

±0.03 

d 2 

-0.665 

±0.005 

-0.5976 

±0.07 

-0.650 

±0.005 

-0.6492 

±0.0071 

0-bar 

17 

±2 

19 

±0.7 

15 

±2 

15.6 

±0.8 

rms 

0.065 


0.065 


0.028 


0.028 



Ganymede, Callisto, and Rhea, respectively. The 3P-HG 
model predicts a forward scattering component to the sur- 
faces of Europa’s leading hemisphere not seen in the 2P- 
HG model solutions. However, for Europa’s trailing hemi- 
sphere the 3P-HG model predicts that the surface is back- 
scattering, while the 2P-HG model predicts an additional 
forward scattering component. Both models show the same 
single scattering function variations with wavelength for 
Ganymede’s trailing hemisphere, where the 0.47-/xm single 
scattering functions are nearly isotropic while the 0.55-jxm 
single scattering functions are highly backscattering. The 
only variations with model seen for Ganymede’s leading 
hemisphere are those for the 0.55-/xm observations. Both 
models show that the single scattering function is backscat- 
tering at 0.47 /xm, but the 2P-HG model predicts that the 
0.55-/xm single scattering function has a forward scattering 
component larger than that predicted by the 3P-HG model. 
However, this may be due to the fact that there were no 
large phase angle observations at 0.55 /xm for the leading 
hemisphere of Ganymede. The 3P-HG model predicts that 
Callisto’s trailing hemisphere is very strongly forward scat- 
tering while the 2P-HG model shows some variations with 
wavelength, where the 0.47-/xm single scattering function 
has both forward and backscattering components and the 
0.55-/xm single scattering function has a low broad back- 
scattering lobe. The 3P-HG model predicts that Callisto’s 
leading hemisphere is nearly isotropic, with a small narrow 
backscattering component, while the 2P-HG model shows 
some variations with wavelength where the 0.47-/xm single 
scattering function is highly forward scattering and has a 
minor backscattering peak and the 0.55-/xm single scatter- 
ing function has a very low, broad backscattering lobe. 
However, the phase curve data for Callisto’s leading hemi- 
sphere have no large phase angle observations, so modeling 
results of this hemisphere need to be interpreted with cau- 
tion. The 2P-HG model predictions for the single scattering 
function of Rhea are similar for both hemispheres at both 
wavelengths. It shows both a forward and backscattering 


component where the forward scattering component is the 
stronger of the two. The 3P-HG model predicts very differ- 
ent scattering functions between the two hemispheres. At 
0.55 gm the predicted single scattering function for the 
leading hemisphere has both components, as seen in the 
2P-HG solutions; however, the forward scattering compo- 
nent is weaker than in the 2P-HG solutions. For the trailing 
hemisphere at 0.55 /xm the 3P-HG function predicts a 
completely backscattering scattering function. At 0.47 /xm 
the 3P-HG function predicts an extremely large forward 
scattered component for Rhea’s leading hemisphere, but a 
moderately forward scattering function for Rhea’s trailing 
hemisphere, similar to that predicted for the 0.55-/xm lead- 
ing hemisphere 3P-HG solution. 

DISCUSSION 

Hapke’s model (Hapke 1981, 1984, 1986) is mathemati- 
cally complex. The changes incurred in a phase curve de- 
rived with Hapke’s model by changing the value of one 
model parameter can often be counterbalanced by chang- 
ing a second parameter. This is why unique solutions to this 
model (especially for disk-integrated observations) require 
data sets that have a broad range of phase angle coverage. 
Large phase angle data (observations greater than 100° 
solar phase) constrain the values of #-bar, w, and the single 
scattering function parameters. Moderate phase angle data 
(observations between 20° and 100° solar phase) combined 
with large phase angle observations also constrain the sin- 
gle particle scattering function, which in turns influences 
the values of 0-bar and w. This is why the results of the 
satellite modeling (of those satellites with adequate larger 
phase angle coverage) that show substantial differences in 
the single scattering function between the 2P-HG Hapke 
model and the 3P-HG Hapke model also show differences 
in w and/or 0-bar between the two model results. In this 
section we focus our discussion on the single particle scat- 
tering results. However, a caveat must be stated. Not all 
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FIG. 8. The single particle scattering functions for the five terrestrial snow data sets where the solid line represents the 2P-HG single scattering 
solution and the dashed line represents the 3P-HG single scattering solution. The graphs where no dashed line is seen indicates that the two single 
scattering functions are equivalent, and the dashed line is overlaid by the solid line. 
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FIG. 9. The single particle scattering functions for Europa’s leading 
and trailing hemispheres, plotted against the Hapke modeling solutions 
using the 2P-HG single scattering function for the 0.47-^m (solid line) 
and 0.55-^tm (short dashed line) observations and the Hapke modeling 
solutions using the 3P-HG single scattering function for the 0.47-yLim (long 
dashed line) and 0.55-,um (dash-dot line) observations. 


the satellite phase curves presented here have the neces- 
sary phase angle coverage to constrain adequately all the 
model parameters. Several are missing large phase angle 
observations, and therefore the values of the single scatter- 
ing parameters along with the surface roughness parameter 
are poorly constrained, as demonstrated by the error bars 
listed in Table V. Domingue and Verbiscer (1997) discuss 
this in greater detail. We have chosen to present these 
phase curves in order to demonstrate the effects of poor 
phase angle coverage on the modeling results. The follow- 
ing discussion is based on the model parameter values 
given in Table II; however, the error bars in Table V 
provide a measure of the accuracy of our interpretations 
based on the available phase angle coverage. 

One of the basic assumptions in Hapke’s model regard- 
ing the scattering behavior of planetary regoliths is that 


while the singly scattered component may be anisotropic, 
the multiply scattered component can be approximated to 
first order as isotropic. As the albedo of a surface increases, 
the role of multiple scattering becomes more important; 
therefore, the assumption of isotropic multiple scattering 
for high albedo surfaces, such as snow and many icy satel- 
lites, is no longer valid (Mishchenko 1994, Goguen 1997, 
Verbiscer and Helfenstein 1997). The model used in this 
study, the results presented in Tables II-V, and the text 
assume that the multiply scattered component is isotropic. 
If the anisotropy in the multiply scattered component is 
not coupled to the anisotropy in the single scattered com- 
ponent, then correcting for anisotropic scattering in the 
multiply scattered component introduces a second scatter- 
ing function and thereby one to three additional parame- 
ters into Hapke’s model. The icy satellite data are not 
extensive enough to constrain the number of parameters 
of such a model. However, if the assumption is made that 
the anisotropy in the single and multiply scattered compo- 
nents can be described by the same scattering function, 
then there are no additional parameters (Verbiscer 1991). 



0 50 100 150 200 

phase angle (degrees) 
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FIG. 10. Same as described in the legend for Fig. 9, for Ganymede. 
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Callisto Leading Hemisphere 



phase angle (degrees) 


Callisto Trailing Hemisphere 



FIG. 11. Same as described in the legend for Fig. 9, for Callisto. 


We refit the snow data using a single-parameter Henyey- 
Greenstein function to describe the anisotropic scattering 
of both the singly and the multiply scattered components. 
Since our initial results of modeling the snow data showed 
that both the 2P-HG and the 3P-HG functions gave solu- 
tions which reduced to a single parameter Henyey- 
Greenstein function we felt that such comparisons would 
be valid. In addition we set bo = 0 since we found that in 
the isotropic multiple scattering case the solutions with 
bo = 0 and bo # 0 showed no substantial differences. 
When we compared our anisotropic multiple scattering 
results to our isotropic multiple scattering results we see 
no difference outside of our error bars for the scattering 
function or the single scattering albedo; however, the rms 
values for the anisotropic case are nearly twice those for 
the isotropic case. Since the anisotropic multiple scattering 
model gives similar results (well within the error bars but 
with higher rms values) to the isotropic model for bright 
terrestrial snows, we hypothesized that the anisotropic 
model would produce similar results as the isotropic model 
for the darker icy satellite surfaces. 


Domingue and Verbiscer (1997) applied both an iso- 
tropic and anisotropic multiple scattering Hapke model to 
the phase curve observations of Europa and Ganymede, 
the brightest of the Galilean satellites. Since Rhea is com- 
parable in albedo to Europa, we also applied an anisotropic 
model (using a 3P-HG scattering function, as was used in 
the Domingue and Verbiscer analysis) to the Rhea leading 
hemisphere phase curves. No application was made to 
Rhea’s trailing hemisphere observations since the large 
phase angle coverage is less than 50°. The anisotropic multi- 
ple scattering modeling results for Rhea are given in Table 
VI and compared to the isotropic multiple scattering re- 
sults. The isotropic versus anisotropic multiple scattering 
modeling results of the icy satellites are commensurate 
with the corresponding modeling results for the terrestrial 
snows. The most substantial differences (differences out- 
side the error bars) between the isotropic and the aniso- 
tropic results for the icy satellites were for phase angle 
ranges not constrained by the data, such as the forward 
scattering direction (a > 120°). Comparisons of the pre- 


Rhea Leading Hemisphere 
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FIG. 13. A plot of the two-dimensional 2P-HG single scattering parameter space ( b 2 versus c 2 ) of Hartman and Domingue (1995) solutions to 
the McGuire and Hapke (1995) artificial particles (■, smooth, clear sphere; #, rough, clear sphere; A, pitted, clear sphere; ♦, irregular shaped 
particle; □, sphere with Iow r density of internal scatterers; O, sphere with moderate density of internal scatterers; A, sphere with high density of 
internal scatterers; O, agglutinate) versus our 2P-HG single scattering parameter results for terrestrial snow (*); Europa (X); Ganymede (0); Callisto 
(ffi), and Rhea (+). 


dieted phase curves and of the single scattering function 
were nearly identical between the isotropic and anisotropic 
model parameters within those phase angles covered by 
the data sets. Any divergences in the predicted phase curve 
behavior or the single scattering function behavior occur 
only at phase angle ranges beyond those constrained by 
observations. There were also significant differences (dif- 
ferences outside of the error bars) in the opposition ampli- 
tude parameter (bo) for both Europa and Ganymede 
(Domingue and Verbiscer 1997). The differences between 
the isotropic bo values and the anisotropic bo values show 
no systematic trend; however, the isotropic values show 
smaller differences in bo between the two wavelength solu- 
tions than between the anisotropic model solutions. The 
parameter values between the isotropic and anisotropic 


multiple scattering models for the icy satellites are qualita- 
tively similar; they predict similar single scattering func- 
tions for phase angle ranges within the observations, simi- 
lar particle opacities (bo < 1 is transparent, bo > 1 is 
opaque), similar surface roughness, and similar single scat- 
tering albedos. The quantitative differences, such as pre- 
dicted single scattering behavior beyond the phase angle 
coverage and the degree of opacity, can be attributed to 
the different mathematical coupling of the scattering func- 
tion to the opposition parameters between the isotropic 
and the anisotropic models. Predictions of scattering be- 
havior for phase angles not contained in a data set are 
going to be strongly model dependent. For the analyses 
presented in this study we have chosen the isotropic multi- 
ple scattering model since the solutions are qualitatively 
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b 3 


FIG. 14. Same as described in the legend for Fig. 13, for the three-dimensional 3P-HG single scattering parameter space defined by b } and c } . 


similar to the anisotropic results and, in addition, the iso- and the satellites. Figures 14 through 16 are similar plots 
tropic solutions have opposition width parameter ( h ) val- of the 3P-HG parameter space. 

ues and macroscopic roughness parameter (0-bar) values Figure 13 shows that our 2P-HG model scattering results 
which are invariant with wavelength, consistent with their suggest that terrestrial snows are most similar in structure 
physical definitions. to the Hartman and Domingue (1997) category A type 

Hartman and Domingue (1997) have shown that there particles. Category A type particles are those which depart 
are no differences in the quality of fit between the 2P-HG the least from a smooth, clear sphere. The frost values, 
and the 3P-HG functions when applied to the McGuire however, plot outside of the laboratory defined areas. Frost 
and Hapke (1995) single scattering laboratory data. Our is shown to have higher values of b 2 than any of the Hart- 
results in modeling the terrestrial snows support their find- man and Domingue particle categories; however, it does 
ings. McGuire and Hapke (1995) and Hartman and Dorn- have similar c 2 values to their category A particles. The 
ingue (1997) plotted the single scattering parameter values 2P-HG values for the icy satellites plot outside of the labo- 
against each other to see if specific particle types cluster ratory defined parameter space. This indicates that there 
in certain areas of single scattering parameter space. The are a variety of particle physical types present on their 
purpose of this was to obtain a mechanism for correlating surfaces and the single scattering parameters represent an 
the single scattering parameters with physical particle average. The 2P-HG values for Europa indicate that the 
structures and textures. Figure 13 is a plot of the 2P-HG structural characteristics of its regolith particles are a corn- 
parameter space that includes the results of Hartman and bination of irregular shape and presence of internal scatter- 
Domingue (1997) with our results for the terrestrial snows ers. These values also indicate a higher presence of internal 
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FIG. 15. Same as described in the legend for Fig. 13, for the three-dimensional 3P-HG single scattering parameter space defined by b$ and d^. 


scatterers for the leading hemisphere with respect to the 
trailing hemisphere. This is consistent with the preferential 
bombardment of the trailing hemisphere by magneto- 
spheric ions (Sack etal 1991, 1992). The 2P-HG results for 
Ganymede and Callisto are scattered and predict different 
particle types dependent on wavelength. The particle types 
predicted should be independent of wavelength. However, 
there was no large phase angle coverage for Ganymede in 
the 0.55-^im observations and none for Callisto's leading 
hemisphere, which may account for some of these inconsis- 
tencies. The Rhea parameters most closely correlate to the 
category A type particles, with possible departures in shape 
from perfect spheres. 

Figures 14 through 16 plot our 3P-HG modeling results 
against the 3P-HG laboratory modeling results of Hartman 
and Domingue (1997). These figures clearly demonstrate 
that the structure of terrestrial snow and the regolith parti- 
cles of the satellites we have chosen to study have nothing 


in common with the structures of the laboratory particles 
of McGuire and Hapke (1995). Our snow results define a 
new area in the parameter space. The parameter values 
for the icy satellites do not fall in this region, supporting 
the Verbiscer and Veverka (1990) hypothesis that the grain 
structures and textures of outer Solar System water ices are 
very different from terrestrial analogs and the laboratory 
particle structures studied by McGuire and Hapke (1995). 
The only laboratory particle which falls into the icy satellite 
region is that of the high albedo (w = 0.67) sphere with a 
high density of internal scatterers. This is a strong indica- 
tion that the particle structures of icy satellite regoliths are 
not uniform, but consist of a variety of particle structure 
and texture types. 

Helfenstein et al (1997) analyzed the McGuire and 
Hapke laboratory data using a 3P-HG function in an effort 
to find an albedo dependence within the particle scattering 
phase functions. In their study Helfenstein et al (1997) 
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FIG. 16. Same as described in the legend for Fig. 13, for the three-dimensional 3P-HG single scattering parameter space defined by d 3 and c 3 . 


define an effective asymmetry parameter, g cff3 , such that 


CONCLUSIONS 


gcfD = <COS 0) = 


/” P 3 (6) cos 6 sin 6 56 
Jlr>(0) sin 6 SO 


= (1 - c 3 )b 3 + c 3 d 3 , 


where 6 is the scattering angle and is related to the phase 
angle by 6 = n - a. For the two parameter HG function 
d y = -b 3 , ^us g c ff2 = b 2 (l ~ 2 c 2 ). Figure 17 plots g cfI3 
for the Hartman and Domingue (1997) laboratory results 
along with our results for the terrestrial snow data and the 
satellite observations. These plots demonstrate that the 
particle structures for terrestrial snow and planetary rego- 
Iiths are much more complicated than the structures of the 
McGuire and Hapke (1995) laboratory particles and that 
the structures of the icy regoliths measured in this study 
are very different than the structure of terrestrial snow 
and frost. 


Our comparisons of the single particle scattering behav- 
ior of terrestrial snows and icy satellite regoliths to the 
laboratory particle scattering measurements of McGuire 
and Hapke (1995) demonstrate that the differences be- 
tween icy satellite regoliths and their terrestrial counter- 
parts are due to particle structures and textures. Terrestrial 
snow particle structures define a region in the single parti- 
cle scattering function parameter space separate from the 
regions defined by the McGuire and Hapke (1995) artificial 
laboratory particles. The particle structures and textures 
of the grains composing icy satellites regoliths are not sim- 
ple or uniform but consist of a variety of particle structure 
and texture types, some of which may be a combination 
of the particle types investigated by McGuire and 
Hapke (1995). 
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FIG. 17 . Plots of the 2P-HG effective scattering parameter (g eff2 ) versus the 3P-HG effective scattering parameter (g eff3 ) for: (top) ETartman 
and Domingue (1995) solutions to the McGuire and Hapke (1995) artificial particles (same symbols as described in the legend to Fig. 13 and 
connected with a solid line) and our satellite modeling results [Europa (X), Ganymede (0), Callisto (ffl), and Rhea (+)]; (bottom) our results for 
terrestrial snow (□, newly fallen snow; A, snow frost; O, rain crusted snow; O, wind blown snow; •, settling snow) versus our satellite modeling 
results (symbols same as described for top graph). 
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